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SOME GENERAL THEOREMS ON THE MOTION OF 
INCOMPRESSIBLE VISCOUS FLUIDS 


GALE YOUNG 


THE UNIVERSITY OF CHICAGO 
(At present: Olivet College, Olivet, Mich.) 


_ Some standard theorems about the motion of single fluids are 
reviewed and extended to the case of several fluids moving through 
each other. Some further results are obtained which do not have a 
counterpart in the case of a single fluid. 


The mechanics of viscous fluids is of importance in the study of 
biological movements, diffusion processes, etc. The purpose of the 
present paper is to collect together some results of a more or less gen- 
eral nature, without entering into any detail as to actual solution of 
the equations of motion. 


Single Fluid 


; 1. An incompressible fluid* can move without viscous dissipa- 
- tion of energy only if it moves as a whole like a rigid body (Lamb, p. 
549). Therefore if it is constrained to have zero velocity over a fixed 
finite area of surface it cannot move at all without dissipating energy 


into heat. 


5 

o . 2. The energy equation for a fluid confined within a fixed bound- 
¢ ary+ at which it has zero velocity is 

f K 

, SRS Ht t SSS Vas (1) 


where V = (u,v,w) is the velocity of the fluid; ¥ = (X,Y,Z) is the 
field force acting on unit volume of fluid; dz denotes the element of 
volume, the integration being throughout the entire region in ques- 
tion; K is the total kinetic energy of the fluid in the region; and f is 
the dissipation. The quantity f is inherently positive, and can vanish 

_ only if the fluid is at rest throughout the region. Equation (1) is per- 

haps sufficiently obvious, but it may be derived by multiplying the 

| * More exactly, a connected mass of such a fluid. Two separate portions of 
fluid can, of course, move relatively to each other without dissipation. 


Here, as throughout, the boundary may consist of one or more internal 
closed surfaces in addition to the external one; i.e. the region it bounds may be 


periphractic (Lamb, p. 38). 
3 145 
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first three equations of (5) below by u, v, w respectively, and then 
adding, integrating, and transforming suitably (Lamb, p. 8). 


If the volume forces have at each instant a potential, X(t) = —V 
Q(t), the last term drops out* leaving simply 
dK 
eee (2) 
dt f 


Thus under the operation of potential forees there is a unique 
steady state, namely that wherein the fluid is at rest throughout the 
region. This is furthermore a stable state, since any imported motion 
dies out to zero by dissipating its kinetic energy into heat. The steady 
state attained is independent of the force field X , provided that at 
each instant it has a potential Q. The dissipation in the steady mo- 
tion is less than in any other motion having the same boundary ve- 
locities. 

3. For a specified motion, K in (2) is proportional to the fluid 
density p, while f is proportional to the viscosity coefficient 7. Thus 
it is seen that the rate at which the motion dies out increases with in- 
creasing u — 7/p. Some idea of this may be gained by supposing that 
the velocity components were to die out uniformly throughout the 
region, i.e. u = U, y(t), etc. Then K and f both vary as y?, and upon 
integrating (2) we obtain 


y= Ce. (3) 


where a > 0 depends only upon the initial velocities and is homogene- 
ous of degree zero in them. Thus u plays somewhat the role of an ex- 
ponential decay factor. It was termed by Maxwell the kinematical 
viscosity. 

4. Next consider the case where non-zero velocities for the fluid 
are prescribed over the boundary of a moving region. These bound- 
ary velocities are not entirely arbitrary, since the fluid motion must 
at each instant satisfy the equation of continuity V - V = 0 through- 
out the region. This requires that the relation 


SfV-ds=0 (4) 


be satisfied identically on the moving boundary. Apart from this re- 
striction the boundary velocities of the fluid may be assigned arbi- 
trarily as functions of time. 


5. The complete equations of motion of an incompressible fluid 
are 


* By transforming to a sur pre integral over ye pounders 
Ti — S, 


and remembering that V — 0 AG the boundary. 
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| where p is the mean pressure of the fluid (Lamb, p. 543), and 
Du. eu 


etc. 

Not much progress can be made with these equations in general. 
For slow motions, however, there are results corresponding to those 
of section 2. If a given motion has its velocities multiplied by « then 
all the terms on the right side of (5) vary as « except those in V- V u 
ete. which vary as «2. Thus as « decreases the latter become of less 
and less importance, and in the limit « = 0, we have the so-called 
equations of slow motion: 


G) 0 
OP Xia Wid disato t 


Ox ot 
oye) Ov 
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6. Let (u,v,w) be one fluid motion, and let (w” =wt+w,v" = 
yt+v',w" =w + w’) be any other motion maintaining the same ve- 
locities on the boundary. Each set of velocity components defines a 
non-negative kinetic energy and a non-negative dissipation. In gen- 
eral, K" # K' + K and f’ # f’' + f. In fact (Korteweg, 1883; Lamb, 
p. 584) the instantaneous dissipations are related by 


fr=ftf—2n SSS WV ut vv + w'Vw) dr. (8) 
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Consider two motions satisfying (7) and taking place under the 
action of force fields X,.and X,. By using (8) twice with first one 
and then the other as the unprimed motion we obtain an equation 


which may be written* as 


SSP + SSS VAX de; (9) 


where V’= V. — V,, X’ = X. — X,, and the left side is the total rate 
of change of K’ within the moving boundary. 

If the difference force field X’ derives at each instant from a po- 
tential this reduces to 


eee (10) 


corresponding to (2). 

Two force fields will be said to be equivalent if, within the mov- 
ing region in question, their difference has at each instant a poten- 
tial. The set of all force fields equivalent to a given one is the equiva- 
lence class of that force field. 

Then from (10) it is seen that any two slow motions maintain- 
ing the same (time variable) velocities over the same (moving) 
boundary, and taking place under the action of equivalent force fields, 
tend to become ultimately identical regardless of different initial con- 
ditions. If at some instant the two motions are identical they will 
remain so thereafter, since K’ can change only by decreasing and it is 
already at its minimum value of zero. It follows that a slow motion 
is uniquely determined by its initial and boundary velocities and the 
equivalence class of its force field. 

7. Equation (10) for the decay of the difference motion is quite 
the same as (2), and the remarks of section 3 apply. Thus the rate 
at which the initial conditions die out in a fluid motion increases with 
increasing uw. If after a certain instant two motions have equivalent 
force fields and the same boundary velocities, they tend to become 
identical at a rate which increases with increasing 4. Thus with in- 
creasing « a motion depends less and less on its past history and be- 


* By an argument involving the fact that V’ is always zero at the moving 
boundary, and that the total time derivative of an integral taken throughout a 
moving region is 


reel ; 


d A lesa 
rie ed lives T+ f/fhW.ds, 


where W is the velocity of the moving boundary. 
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comes more and more nearly determined by the instantaneous bound- 
ary velocities and force class. 

The limit of this trend is a motion which at each instant satisfies 
the equations obtained from (5) by omitting the terms in p, namely 


é 


6) as 

0p 

eS § ei 2 

oe, EN ae 
(11) 

op 
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For two motions satisfying these equations and having equivalent 
force fields and the same boundary velocities the use of (8) twice 
gives simply f’ = 0, so that the motions must be identical. It follows 
that a solution of (11) is uniquely determined by the boundary ve- 
locities and the class of the force field. 

Let V be a motion satisfying (7), and let V. = V —V' be a mo- 
tion with the same boundary velocities, and equivalent force field, and 
at each instant satisfying (11). Using (8) twice then gives 

dK’ ; , OVo 
Car =—f Pp af shy f V'- at 

If V satisfies (5) rather than (7) there are additional terms hav- 
ing a factor p. 

As p/7 goes to zero so does f'/y in (12). Thus with increasing u 
a motion tends to becomes identical with Vo, which depends only up- 
on the instantaneous boundary velocities and force class. 

8 Jt is obvious that there can be no steady motion in a fixed 
region unless the boundary velocities are maintained steady in time. 
From (5) it is seen that a further necessary condition for steady mo- 
tion is that X — dp/ox, etc., be independent of time; i.e. that X be 
equivalent to a steady force field. 

But if these conditions are satisfied Vo as defined above is steady 


and (12) reduces to 


Gir. (12) 


dk’ 
a i | 13 
di 1) (13) 


so that a slow motion with assigned steady boundary velocities and 
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steady force class tends stably to a unique steady state. This follows 
otherwise from considerations given in section 10 below. 

9. Let V, be a slow, steady, potential-force motion and V be 
any other motion with the same boundary velocities. Then (8) gives 


f=f+f, (14) 


where V’'= V — V,. Thus the dissipation in slow, steady motion un- 
der the operation of forces having at each instant a potential is less 
than in any other motion having the same boundary velocities. 

Some special motions satisfying the complete equations of mo- 
tion (5) are known (Rayleigh, 1913; Lamb, p. 585) which also make 
the dissipation a minimum; in general, however, any actual motion 
will have a higher dissipation than the corresponding potential force 
solution of (11). 

More generally (8) and (11) yield 


PALS falc Kade =F 25S Lo Vek dette, (15) 


so that the excess of the dissipation over twice the rate at which the 
volume forces are doing work upon the fluid is less in slow steady mo- 
tion than in any other motion having the same boundary velocities 
and the same force field. 

10. From sections 8 and 9 it is seen that with steady boundary 
velocities and potential forces a slow motion tends to a steady state 
of minimum dissipation. It can be further shown that the dissipation 
decreases monotonically in this process. 

For any motion with steady boundary velocities we have (Lamb, 
p. 585) 

df , : : 

tal —2nSfSff (uV2u + vV20 + wV2w) dz, (16) 
where the dots denote partial time derivatives. If the motion satisfies 
(7) with potential forces, this becomes 


—=—2pSSf (w+ ot w) de, (17) 
so that f continually decreases until the motion is steady. 
More generally for a slow motion under forces equivalent to a 


steady field X the left side of (17) is replaced by be (f — 29), where 


9=SSfV-Xd. (18) 


This describes the monotonic decrease of the quantity which was 
shown in (15) to attain its minimum value in the steady state. 
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Several Fluids 


11. In connection with diffusion processes it is appropriate 
(Young, 1938) to consider the motion of different fluids through each 
other, on the supposition that they exert frictional drag forces on 
each other whenever there is relative motion between them. Thus re- 
garded, Fick’s law describes the diffusion of a solute as the motion 
of an incompressible fluid whose inertial and viscous forces are neg- 
lected in comparison with the drag force exerted by the solvent. 

To the extent to which the diffusion coefficient is constant, Fick’s 
law makes the drag force between solvent and solute proportional to 
their relative velocities. This is an assumption that in one form or 
another has had considerable acceptance and support (Fletcher, 1911; 
Smoluchowski, 1916; Burger, 1918; Ehrenfest, 1918, Wiener, 1921; 
Weyssenhoff, 1925; Chapman, 1928; Young, 1938), and it will be 
made in what follows here. We suppose, therefore, that the volume 
drag force exerted by fluid 7 on fluid 7 is given by 


Fi,=hij (Vi — Vi), (19) 


where the k;; are positive constants. For the mutual drag forces to 
be equal and opposite it is necessary that ki; = k;,;. Note that when 
a solute is pictured as an incompressible fluid the velocity is propor- 
tional to the mass rate of flow, i.e. to what would ordinarily be de- 
noted by the product of the density and the mean drift velocity of the 
solute particles. The total volume force on fluid 7 is 


iat AG (20) 
j 


where X; is the resultant of all other forces besides the drag forces. 
For convenience X; will be referred to as the external force on fluid 7. 

12. Using the total force (20) in equation (1) and summing 
for all the fluids gives the total energy equation for a region in which 
all fluids have zero boundary velocities; namely 


SSK HFS SSSXVide, (21) 
where 
F= Sf t+ Dd Oy 
oe (22) 
O,,=4 SSS Fis: Vi —Vi) de. 
F is the total dissipation within the system due to viscous and fric- 


tional drag forces. It is inherently positive, as is apparent from (19) : 
and vanishes only when all the fluids are at rest. For this conclusion 
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it would have been sufficient to suppose merely that the scalar prod- 
uct under the integral in ®;; is non-negative, instead of the more re- 
strictive assumption (19). This, however, would not suffice for sec- 
tions 13 and 15 below which assume (19) specifically. 

If the external forces X; have at each instant a potential , (21) 


reduces to 
d 
Sk. = re. 23 
dt ee 


This is analagous to (2), and shows that a system with zero bound- 
ary velocities and potential external force fields tends stably to a 
unique steady state of zero motion and minimum total dissipation F’. 

Writing out (5) for each component fluid shows that steady zero 
motion requires X; = 0p;/dx, etc., so that the system can be at rest 
only if all the external force fields are equivalent to zero. 

Thus the slow motion of one or more incompressible fluids with 
relative drag forces as given by (19) is uniquely determined by the 
initial and boundary velocities and the equivalence classes of the ex- 
ternal force fields. Two motions with different initial velocities but 
with identical boundary velocities and equivalent external force fields 
tend to become ultimately identical. If the boundary velocities and 
the equivalence class of the external force field are maintained steady 
in time then the motion tends stably to a unique steady state. 

13. Using the total force as given by (20) and (19) in (9) 
gives for the difference of two slow motions maintaining the same 
boundary velocities 


d 
Bueis P+ SSffxiVin, (24) 


so that the conclusions of section 6 apply to a system of fluids moving 
through each other under the influence of drag forces (19), just as to 
the motion of a single fluid. 


For slow motion with steady boundary velocities and potential ex- 
ternal fields (16) gives 


dF ' , : 
—=—2a a SSS ut foe + wer) ar: (25) 
dt 
or if the external fields are equivalent to steady fields, which is neces- 
sary for existence of a steady motion, the left side is replaced by 


d : 
ag {F — 2G) where G is the sum of g; in (18). Thus the conclusion of 


section 10 also extends to the case of multiple fluids with drag forces 
as in (19). 
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14. The total stress system due to several moving fluids is the 
sum of the stress systems for each fluid, which in turn involve sums 
of the pressures p; and space derivatives of ; times the various 
velocity components (Lamb, p. 544). The total stress thus depends 
only upon the quantities 


P=2Dpi 


y — poyP Vi; . 


(26) 


Writing the equations of slow steady motion (11) for each fluid, using 
(20), and summing over the fluids gives 


oP 


oe ae Le 
Oe eis, 
oy 
(27) 
SEO ieee 
0z 


OU ov ow 
Se ry ee | 
where P and (u, v, w) = V are the quantities defined in (26), and X 
is the sum of the external force fields X; of (20). But (27) describes 
the slow steady state motion of a single fluid under a force field X, 
which is determined by the boundary velocities. Hence the total stress 
ean be found without determining the individual fluid velocities V; . 
In particular if X is equivalent to zero then the total stress is merely 
the sum of the stresses of each fluid in steady force-free motion in 
the absence of the other fluids and with its own boundary velocities. 

If all the fluids have the same kinematical viscosity mu a corres- 
ponding result may be obtained from the time variable slow motion 
equations (7). This results in the addition of terms — (1/1) (0u/ot), 
etc., on the right sides of the first three equations in (27), and reduces 
the calculation of the total stress in a varying multiple fluid system 
to that of a single fluid motion. | 

Note that the results of this article involve only the assumption 
that F,;—= — Fj: ; they do not depend upon (19). 

15. Now suppose that one of the fluids (say No. 1) in a system 
is identically at rest. Assuming the specialized drag forces (19) and 
introducing the total force (20) into the complete equations of motion 
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(5) for a single fluid we see that this implies the equations 


suas X, = Dj bay Uy 


Ox 

OP. _ y= 3 deus 0; (28) 
oy 

CPs 21 = Dj bj W; , 

0zZ 


where X,= (X,, Y:, Z;) is the external force acting on the station- 
ary fluid. Differentiating and adding and remembering that V -V; 
= 0 for each fluid gives 


V*pi=V-X. (29) 
If X, has at each instant a potential (29) becomes 
V?P=0, (30) 
where 
P= pao, 
(31) 
AX,=- VQ, > 
while (28) becomes 
where 
V =; 4; V;. (33) 


But then P is determined by the value of its normal derivative over 
the boundary, i.e., by the normal components of the fluid velocities. 
The tangential velocities can not be assigned arbitrarily in addition 
without rendering (28) inconsistent. Hence a slight change in the 
boundary velocities of any of the other fluids can enforce motion on 
the part of fluid No. 1, and except in special cases the motion of any 
fluid entails that of all the fluids. A special case of this result has 
been noted previously (Young, 1940). 

If X, does not have a potential it is still seen from (29) that D1 
is determined to within a harmonic function; 


D=at+H, (34) 


where a is determined by X, and H is an arbitrary harmonic. Then 
(28) becomes 


eA Vas (35) 
and the argument goes through as before. 


a eS ee 
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A NEURAL MECHANISM FOR DISCRIMINATION II: 
VISUALLY PERCEIVED LENGTHS AND DISTANCES 


ALSTON S. HOUSEHOLDER 
THE UNIVERSITY OF CHICAGO 


; A previously discussed neural mechanism for the discrimi- 
nation of intensities is here applied to the judgment of visual lengths 
and distances on the assumption that the “intensity” associated with 
the magnitude being perceived is the intensity of innervation of the 
appropriate eye muscles necessary for scanning and fixating. Com- 


parison with experimental data is made in the case of the judgment 
of lengths. 


1. The problem. In previous papers (Householder 1939, 1940: 
these will be denoted hereafter by DI and DII.) a hypothetical neural 
mechanism was described which has the property of diverting to 
different neural pathways the nervous responses from stimuli of dif- 
ferent intensities. Furthermore, from this mechanism as it was set 
up it was possible to deduce the form of the Weber function that 
would result (i.e. the function describing the stimulus which is “just 
noticeably greater” than a given stimulus) and so to compare this 
theoretical function with those obtained experimentally. Comparisons 
were made with empirical data taken from the literature for visual 
brightness, for loudness, for tactile pressures, and for lifted weights. 
For the first three types of data the intensity of the stimulus S acting 
upon the mechanism proper was supposed proportional to the inten- 
sity of the peripheral stimulus whereas for lifted weights account was 
taken of the probable effect of the muscle receptors in transmitting 
the stimulus. 

In the present investigation we wish to consider how the same 
type of mechanism might operate in the making of judgments on spa- 
tial magnitudes visually perceived. Inasmuch as the mechanism as- 
sumes that stimuli of graded intensity are applied, it is necessary to 
tell what physiological intensities are supposed to yield the experi- 
ence of the extension of the stimulus-object, for only if this is pos- 
sible will the mechanism in question be applicable, and the derivation 
of the associated Weber function be legitimate. Before doing this, 
however, it is convenient to introduce an approximation to the basic 
formula derived in DI, which gives the Weber function in terms of 
the intensity acting upon the mechanism itself. 
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2. An approximation. If u is proportional to the intensity of 
the stimulus acting upon the mechanism and if wu and u(1+6) are the 
intensities corresponding to “just discriminably different” peripheral 


stimuli, then according to equations (16), (22), and (23) from DI 
uz? — 2(a+1)2— (#+1)?=0 
ux? + x? -1=0 (1) 
x+2z=—6 


where x and z are auxiliary variables. The solution giving 6 as a 
function of u can actually be obtained in closed form; this was done 
in DI, and it was further shown that 6 is monotonically decreasing 
from o for u=0 to 0 for w= o. Nevertheless, this closed solution 
is too complex and unmanageable. It is simpler to obtain uw and 6 as 
functions of a single parameter, and from this parametric represen- 
tation a much simpler but still adequate approximate expression is 
readily derived. 
From the second of equations (1) we have 


_ eae 
LL ene 


(2) 


Hence if we set this expression for uw and substitute 
2— ry 


into the first of these equations we obtain an equation which readily 
simplifies to the form 


Ql 0) Ary ee (xe 1) — Oe 


after suppressing a factor «+1 which is necessarily positive. From 
this we have 


ta ue 
yy +2y+1° (3) 
Evidently 
d=a(1+y), 
or 


MG ral) 49 e 
y+ aytl (4) 


Thus if we substitute the expression (3) for x into (2) we have 


wa t¥ Yt DY +3) + 2y +1) 


(y3 —_ 1)3 ? (5) 


zr TC 


ee em ha 
. 
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and (4) and (5) together constitute a parametric representation as 
desired. 

However we note that 6 is monotonically increasing and-w mo- 
notonically decreasing in y , and that the range of y for which wu and 
6 are positive, and hence meaningful, is from y= 1 to y= 0. (The 
proof of this statement is contained in the ensuing discussion.) It is 
more convenient therefore, to introduce the new parameter f, where 

= ee ; (6) 
so that wu is monotonically increasing, 6 monotonically decreasing in 
t, and the range of tis fromt—0Otot=@. This gives 


pA +1) (2t + 1) (202 + 2t + 1) (40 + 5 + 36+ 1) 
(8S ot 1)= 


? 


(7) 
OAR UN Ch ore 3t4+1) 
— ¢(4@+524+3t4+1) © 
Let us make the substitutions 
n—logo, E=logu. (8) 


Then we find 
dyn/dt = —(24t + T20 + Q7tt + 54t? + 210? + 6 +1)/ 
fé(2t + 1) (8@ + 8641) 46 + 5t2+3¢4+1)], 
dé/dt = (144° + 648¢ + 13500" + 1704t® + 1443¢* + 858% 
4+ 36323 + 108¢2 + 21¢ + 2)/ [t(t + 1) (2t + MGA Sees) 
x (gt? + 3t +1) (40 + 52 + 38f+1)]. 


The monotonicity of € and 7, and hence of wu and 6 as asserted above, 
is evident from these two formulas. From these we find by dividing 
that 


dn/dé = — (48t° + 2408 + 5347 + 696t® + 591t® + 345¢ + 143% 
+ 43t2 + 9t + 1)/(144t + 64828 + 135007 + 1704t°® (9) 
+ 1443¢ + 858t¢ + 363¢* + 1082 + 21t + 2) 


=—1/(2+t7,)=—(1+T:)/8, 


where 
7 == (AS + 1680 +- 2g0t6 + 312t + 2614 + 168¢% + 77t2 + 22t + 3)/ 


(4gt? + 2400 + 534t7 + 6960 + 591t + 345t4 
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+ 143 + 43¢ + 9¢+ 1) 


and 
T= (T2t® + 252t7 + 384¢* + 3300 + 177t* + 660 + 21t? +.6¢( + 1)/ 


(144° + 648¢8 + 1350¢7 + 1704¢* + 1443¢° + 858¢* 
+ 36323 + 108¢? + 21¢ + 2). 


Since T, and T> are always positive when ¢ is positive, it is evident 
that 


1/2 Sdy/dé = 1/8, (10) 


and that this derivative takes the extreme values —1/2 and —1/3 only 
when t = 0 and when t = %, respectively. Because of the slight vari- 
ation in the value of this derivative, 1 can be represented with good 
approximation over a finite range as a linear function of € , and hence 
6 can be represented as proportional to a power of wu over this same 
range. The representation can be determined by taking the first term 
in the Taylor series expansion of 7 — 7) in powers of € — &, for some 
& and y., or, perhaps better, by taking the linear interpolation for- 
mula determined by some two points (&, 70) and (&,, 7:). The cal- 
culations can be made most conveniently from the simpler expressions 
(4) and (5) in terms of y. 

In any case, we have, for the approximate representation of 6 
as a function of u the expression 


= p4e* (11) 


where f and a are to be determined according to the range over which 
the approximation is required, and where 
Accel 
3= a= >" (12) 
At this point we may point out the similarity of the form (11) 
to the psychophysical law of Cattell and the generalizations thereof 
(Woodworth, 1914; Guilford,1932). It is true that w is not supposed 
to be the intensity of the peripheral stimulus, but to be the amount of 
the excitatory “factor”, i.e. the intensity of the stimulus, which is act- 
ing upon the central discrimination mechanism itself. Nevertheless 
it may not be too unreasonable to suppose that this central intensity 
is roughly proportional to the peripheral intensity for at least some 
sense modalities, in which case equation (11) can be used directly 
with wu as the peripheral intensity. This was done in the paper DI for 
visual brightness, for loudness, and for pressure. 
It should be noted further that a and f are not independent para- 
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meters, and that the only variable parameter actually involved in for- 
mula (11) is the factor of proportionality contained in the variable w . 

3. Visual lengths. It has been remarked above that if the theory 
of intensity discrimination is to be applicable to the discrimination of 
visually perceived magnitudes, then we must suppose the experience 
of extension to be due to some stimulus whose gradations in intensity 
correspond to the magnitude of the extension. This physiological in- 
tensity of stimulation must be supposed to underly and to result in 
the psychological experience of extension in space. 

This physiological intensity we shall assume to be associated with 
the extrinsic eye-muscles. In view of the known distribution of the 
rods and cones on the retina (Hecht, 1938), it does not seem reason- 
able that the required intensity should be due to a summation of the 
responses resulting from stimulation of these elements, though this 
would be otherwise the most natural assumption to make. The ques- 
tion of whether the experience of spatial extension in visual percep- 
tion can be attributed to the role of the eye-muscles has been much 
discussed in qualitative and philosophical terms, and we do not need 
to go into this here (see, e.g., Baird, 1903; Roelofs, 1935, among 
others). In this connection we might mention however, some work of 
Kennard’s (1939) in which lesions in the motor cortex in monkeys 
led to visual disturbances. It is of interest to note, also, that Mach 
(1861) formulated a quantitative theory for the discrimination of 
visually perceived lengths and angles based on the assumption that 
the experience of extension derives from the action of the eye-muscles. 

Recently Rashevsky (1938b, 1940) developed a theory of the per- 
ception of visual forms based on such an assumption. As he points 
out, it cannot be the actual movements of the eye that are responsible 
for the experience of the extension, or for the perception of the form, 
since it is known that the eye does not in fact, or at least need not, 
accurately scan the outlines of the form being perceived. But we may 
assume that the actual movements of the eyes are the resultants of 
innervations from various centers, one of which, if the others were 
inactive, would cause the eye actually to scan the object. This is the 
center Cs in Figure 1, taken from the above mentioned papers. 

Suppose, now, that there is presented for visual perception a sin- 
gle straight-line segment, placed in the horizontal plane of the two 
eyes. For simplicity we consider only one eye. Suppose the segment 
is so placed that one endpoint can be fixated when the eye is in the 
primary position, the segment and the line of sight being then at right 
angles. The innervation from the center Cs is then supposed to be 
such that it tends to cause the eye to rotate so as to fixate the other 
endpoint of the segment. The center Cy, is supposed to be connected 
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with a discriminating mechanism of the type discussed in DI, so that 
the intensity S, or u delivered to this mechanism is proportional to 
the intensity of the innervation from Cs. 

Since the segment is supposed to be horizontal, the rotation of 
the eye, if Cs alone were acting, would be in the horizontal plane, and 
only the action of the horizontal rectus muscles need be considered. 
For vertical, and certainly for oblique, rotations the situation is some- 
what more complicated since more than one pair of muscles is in- 
volved. In the present case the linear contraction of the appropriate, 
lateral or medial, rectus muscle is proportional to the angular rota- 
tion of the bulbus. It may be that the relative, rather than the linear, 
contraction of the muscle should be considered, but the difference 
could amount, at most, to a change of a few per cent in the Weber 
ratio for longer segments. 

If S is the amount of « — 7 produced at the synapse as a result of 
the activity of Cs, and if S is the amount produced as a result of the 


activity of Cy, then S + S is the total amount acting on the fibers 
leading to the muscles, and S is the amount acting on the fibers lead- 
ing to the discrimination mechanism. The relation between the ac- 
tual contraction of the muscle fibers and the quantity S, when S=0, 
depends, in part, upon the number of stages in the sequence of fibers 
between this and the muscle fibers themselves; likewise the relation 
between S and the measure u of « — 7 acting upon the discrimination 
mechanism depends upon the number of stages in this intervening 
sequence. If FE represents the intensity of excitation of a nerve fiber 
(see Rashevsky 1938a or 1940 for a definition of this term), or if it 
represents the linear contraction of the muscle, and S the amount of 
e — 7 acting upon it, we assume that in either case E is the same func- 
tion of S , and hence if the same number of synapses intervene in the 
sequence between Cy and the muscle as between Cy, and the mechan- 
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ism, then we may suppose that the actual rotations of the bulbus is 
given by 


6=a(S+S), (13) 


with, now, S proportional to wu. 

_ If the centers Cs and Cy are completely independent, then for 
scanning a segment which subtends an angle 6, S should be given by 
(13) with S = 0. But if these are not completely independent, S 
should be given by (13) with S taken as some non-zero, positive or 
negative quantity. This may be thought of as corresponding to a cer- 
tain average contribution of the center Cy. No suggestion is made 
here as to the precise manner in which such an interconnection might 
oecur, but certainly the fact that different centers mutually affect one 
ahother is too well known to require discussion. Neither do we sug- 
gest a mechanism for the center Cs itself, though it seems reasonable 
to suppose that its operation depends in some fashion upon condi- 
tioning. 

Instead of (13), various other forms might be suggested and 
would be a priori possible. If the sequence to the muscles is the longer 
sequence, @ should be a function of w which is concave downward; if 
the reverse it should be concave upward. The relation (13) may be 
regarded as a first approximation to either one, and in the absence of 
anatomical evidence, one should be able to decide between the two 
cases by an adequate series of discrimination measurements. 

If we assume (13), then, we have 


u=y(6 — 4) ’ (14) 


by an obvious change of parameter, since u is proportional to S. 
Hence from (11) 


ut+ du= u + pur-*= (6 + 46 — 4) , 


where A is the Weber ratio measured in terms of the angles subtend- 
ed by the discriminated lengths. Hence, after making a further change 


of parameter, 
a Raa d (eet Facial (15) 


In this equation there are two parameters, x and 4). The term o. 
in the exponent is not a parameter, as was noted above, but has the 
value that enables us to approximate to the function 6(u), over the 
range covered, by an expression of the form (11). 

This equation (15) has been compared with the experimental re- 
sults of Chodin (1877) as shown in Figure 2, where the points repre- 
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sent experimental values, and the curve is the graph of equation (15). 
Chodin laid the segments to be compared end-to-end (whence our as- 
sumption that one end-point was to be fixated by the eye in the pri- 
mary position), and obtained J.N.D.’s for lengths of 2.5, 5, 10, 20, 
40, 80 and 160 mm. at a distance of 350 mm. He ran two series of 
experiments, and we have taken, for our comparison, the means of 
these two sets of results. The values of the parameters and of a are 
given on the graph. 

For very short lengths the mechanism necessarily breaks down, 
since according to (15), A would have to be zero for 6 = 6, which 
is obviously absurd. It is not to be supposed that objects whose ret- 
inal images are of the order of the diameter of the fovea should be 
judged by means of the same mechanism. Quite possibly a summa- 
tion of retinal stimulation may be involved for such magnitudes. 

Moreover, for longer segments the assumption that the horizon- 
tal rectus muscles are alone acting becomes less accurate, and devia- 
tions are to be expected here. Whether their activity should “confuse” 
the judgment and make the ratio larger, or should “sharpen” it, 
would be a question to be decided empirically, and the mechanism | 
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would require further elaboration accordingly. It may, indeed, be 
noted that if one attempts to fit only the first four points (for the 
four shorter lengths), the parameters can be chosen to give a very 
accurate representation, in which case, however, the empirical Weber 
ratio for the three longer segments are considerably greater than the 
predicted values. This seems to indicate some kind of a “confusion”, 


rendering the discrimination less fine. But if one supposes that S is 
an increasing function of the innervation to these other muscles then 
this is precisely what one should expect. 

4. Binocular distance judgments. In the judgment of distance 
(from subject to object) in a given situation, obviously a variety of 
cues is possible and any one or any combination of these might ac- 
tually be employed. However, one of these must be due to the degree 
of convergence of the eyes, or to the disparity of the retinal images of 
the object, when this is seen with both eyes. Evidently accommoda- 
tion also plays a part, and the relative importance of these two fac- 
tors might well vary from individual to individual (Carr and Allen, 
1906) and also from situation to situation. 

Suppose that, by some means, accommodation cues, size cues, 
shadings and the like can be eliminated, and the subject is required to 
compare the distances of two objects, successively presented, when 
no cues are possible except the degrees of convergence. We shall sup- 
pose that fusion of the images is accomplished in every case, though 
this assumption is not necessary if we suppose that, as with the pre- 
vious mechanism, the innervation from a certain Cy, produces the w, 
rather than the actual fusional movements themselves. Presumably 
u should be zero when the eyes are in the primary position, and should 
increase as the angle between the visual lines for the two eyes in- 
creases. The simplest assumption would be that u is simply propor- 
tional to this angle 26. On the other hand, however, it is impossible 
to converge the eyes beyond a certain maximal angle, whence a more 
plausible assumption would be that u becomes infinite for some value 
of @ which is less than, or at most equal to 2/2. The simplest function 
of this type is 


u=y 6/ (A — 0). (16) 


_ However, if we suppose that 6, = 2/2 and that convergence is possible 
only up a certain maximal value of uw, then we may write 


u=y tan dé (17) 


as a function very closely approximating to (16). This, however, does 
not deviate appreciably from the linearity assumption, since the graph 
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of the tangent can be fairly closely approximated by a straight line 
for angles up to about 2/4, which is close to the maximal convergence 
possible. 

If we consider (17), and if 2a is the distance between the eyes, 
and x the distance to the object, then evidently 


a/x—tand=u/y, (18) 


whence wu is inversely proportional to the distance. In this case we 
may take the inverse of the distance as the measure of the physiologi- 
cal intensity, and the theoretical Weber ratio so measured is given 
by the function 6(w) . The only parameter involved is the parameter 
ya. 

If the form (16) is assumed, A can be calculated as a function of 
6 from the approximation (11), but it is perhaps simpler to calculate 
the empirical u/y and (1 + 6) u/y from (16), from which 6 can be 
obtained and compared with the theoretical curve. Here an additional 
parameter 6, is introduced. Experiments are now being conducted 
for the purpose of obtaining a series of empirical measurements for 
testing the above formulation. 

This work was aided by a grant from the Dr. Wallace C. and 
Clara A. Abbott Memorial Fund of the University of Chicago. 
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DEFORMATION OF SHELL-SHAPED CELLULAR AGGRE- 
GATES: APPLICATION TO GASTRULATION 


N. RASHEVSKY 
THE UNIVERSITY OF CHICAGO 


: An equation derived in a previous paper is integrated, and ap- 
plied to the theory of invagination of a gastrula. Comparison with 
experimental data is made. 


In a previous paper (Rashevsky, 1940 b) we have derived the 
following equation for the deformation of a cellular shell: 
1 pr 
we dE auhbn 3 ou 
The notations used here are the same as in loc. cit. We shall first 
study the above equation for the case of a constant p , assuming that 
during the process of closure the shell remains spherical. Such a case 
will arise when the cells of the aggregate have a constant polarity 
(Rashevsky, 1940 a). 
Referring to Figure 1, we have again, as in loc cit., the relations: 


FIGURE 1 


V=2ar2(1—cos¢)d; rr =rsing. (2) 


Expressing (1) the quantities 7 and 7’ in terms of ¢ and V, and put- 


ting 
V « 
fu) 3 
P Jaa at a 
16 


9 
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we obtain from equation (1): 


op ee ae (4) 
dt (1—c0s'6) 2? 


Introducing 
—e? x 
en (5) 
and separating the variables, we find 
sin? ade _ pay B=A/23? . (6) 
COs & 
Integrating, we find: 
1 belassine 
Se te lope Bee (7) 
2 sing 


where C is an integration constant. 

The next case which we investigate is that of a polarization vary- 
ing in the manner discussed previously (Rashevsky, 1940 b). To this 
end we should introduce equation (8) of loc. cit. into equation (10) 
of loc. cit. and then the latter into equation (11). The expression for 
p thus obtained should then be introduced into equation (1) of the 
present paper. Such a procedure however leads to very clumsy ex- 
pressions which cannot be handled analytically. We therefore use a 
different approximation method from the one used in loc. cit. in order 
to obtain the desired result. 

Consider first a completely closed spherical shell, whose thick- 
ness 6 is small as compared with its radius 7. We may calculate the 
difference c, — c, in that case by remembering that the amount 42726q 
of substance consumed in the shell per second is equal to the amount 
flowing in. The average amount flowing radially through the shell is 
equal to 4n77(c. — ¢,)Di/d. Since at the inner surface of the shell 
the radial flow is zero, therefore the maximum flow inward at the 
outer surface is 827(c, — ¢,)Di/d. Equating this amount to that 
consumed, we find: 


a= G ere g (8) 
2D; 

An expression of the same form, except for a difference of the 
numerical factor is obtained by introducing equation (8) of loc. cit. 
into equation (10) of loc. cit. after putting in both of them s, = 0 : 
S, = 2, which corresponds to a closed shell. 
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Consider again a completely closed shell, in which, however, one 
part is characterized by a rate of consumption q, , and a diffusion co- 
efficient D;, , while the other is characterized by corresponding values 
gq. and D;.. Since we are interested in the average value of © —¢, 
we may calculate that average value by using equation (8), in which 
we put for q and D; some average values q and D;. 

An open shell, such as shown in cross section on Figure 1 may be 
considered as a closed one, in which a part 2a7?(1 — cos ¢) is char- 
acterized by the values g and D; , while the other part is characterized 
by values 0 and D.. Then we have q= q(1 — cos ¢) /2 and 


D;= [Di(1 — cos ¢) + D-(1 + cos $)]/2. 


Introducing those expressions into equation (8) of this paper 
and the latter into equation (11) of loc. cit., we shall find that p now 
contains a factor 


1—cos¢ 


Di = cos $) + D.(d + 00s 4) (9) 


Instead of equation (4) we now have an equation of the form: 


dp _ A, sin ¢ (10) 
dt  /1—cos¢[Di(1— cos ¢) + D.(1 + cos $)] 


Equation (10) gives, when integrated, 


HrSpies 1+ si 
Ee op ee Bite le (11) 
D 2, 1—sinz 


where again « = ¢/2. 

Such equations may describe approximately the closure of a half 
blastula, as has been suggested (Rashevsky, 1940 a). A similar mech- 
anism may be responsible for the invagination of a gastrula. In the 
latter case things are complicated by the circumstance that the in- 
vagination of the entodermal cell layer is necessarily accompanied by 
a deformation of the ectodermal layer, and the simple equation (1) 
will therefore not hold. If we assume however that at the time of gas- 
trulation the actual polarization forces in the ectoderm are much 
smaller than in the entoderm, then approximately, the effect of the 
ectoderm, which in this case is merely passively “dragged”, will be to 
decrease the right hand side of (1), or in other words to reduce B or B; . 

It may also be questioned as to whether considerations about “‘in- 
duced polarization”, which lead to equation (11) can be applied at 
all to an invaginating gastrula. A somewhat similar situation may 
however arise, when a substance produced by the ectoderm is con- 
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FIGURE 2 


sumed by the entoderm, resulting in a flow as shown schematically 
on Figure 2. The concentration at a would be less than at b , produc- 
ing a polarization which is a function of ¢. Such an assumption does 
not seem implausible at all in view of J. Brachet’s (1934) findings, 
that gastrulation can go on in absence of oxygen, which is the only 
substance with negative metabolism, used in an early embryo in 
quantities sufficient to produce appreciable mechanical forces. A 
transport of a substance within the embryo, as shown on Figure 2, 
would account for such a situation. 

On the closure of half blastulae no quantitative data are yet 
available. Very meager data are available on gastrulation. We shall 
apply equations (7) and (11) to those data, mainly to dllustrate a 
possible quantitative approach to embryological problems. 

One characteristic difficulty arises in this connection. Whereas 
the equations are derived under the assumption that the cellular shell 
remains spherical during the process of deformation, sections of in- 
vaginating gastrula show that those deviate very much from a sphere. 
This difficulty may be overcome in the following way. 

Denoting the total length of the inner periphery of the meridional 
cross section by s — in other words, putting s=2A0’' in Figure 1 — 
we have 


s=2¢r. (12) 
Putting | = 27”, we find from equation (12) and the second equation 


(2): 
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(13) 


If the shell, represented on Figure 1 is deformed without stretching 
or volume change, then it may become rather irregular in shape, and 
the notion of the angle ¢ will lose its meaning. However the ratio I/s 
will remain unchanged, as long as the opening of the shell remains 
approximately circular. 

Considering therefore an actual cellular shell, whose opening is 
approximately circular, but whose meridional cross section is not cir- 
cular, we may define a quantity ¢ as the root of the transcendental 
equation (13). To the extent that our simplifying assumptions hold, 
we may expect this quantity ¢ to follow approximately either equation 
try or CLI). 

On Figure 3 are shown tracings of meridional cross sections of six 
different stages of Amphioxus, taken from a paper E. G. Conklin 
(1932, plate 13). The inner perimeters of the cross section of the 
entoderm, which are shown on Figure 2 in somewhat heavier lines, 
were measured with a map measure. The dimensions of Figure 2 are 
considerably reduced, compared with the original. The wheel of the 
map measure was run over the same perimeter 10 times, without re- 
setting to zero, and the final reading divided by 10. The distances be- 
tween the ends of the heavier lines were taken as values of |. The com- 
puted values of ¢ for each stage are given in radians underneath each 
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FIGURE 4 


section. Above each section is given the time in hours, taken from 
Conklin’s paper. 

The simplest check of our equations would be to plot the left hand 
sides of equation (7) or (11) against ¢ ; this should give us a straight 
line. This method however is too sensitive for such inaccurate data 
as we have. In the neighborhood of ¢ = 12, that is for the later stages, 
an error of 8% in the value of ¢ results in an error of 100% in the 
value of the left hand sides of (7) or (11). 

Therefore experimental values of t were plotted directly against 
¢, and the plots turned by 90°, so as to give ¢ against t. The results 
are shown on Figure 4, the dotted line representing equation (7), 
while the full line represents equation (11). The agreement of the 
latter with the experimental values is fair. The value of (D.— D;)/Di, 
was taken as 1.04. 

It must be once more emphasized, that no definite conclusion can 
be drawn yet from such a comparison, in view of the paucity of the 
data. But such considerations are very important in principle, since 
they show a way to new quantitative measurements in embryology. 

If instead of five different stages, we had fifteen, and if for each 
stage we had a large number of cross sections, so that accurate aver- 
age values of ¢ could be calculated, a comparison of theory and ex- 
periment would become very meaningful. It must be remarked that 
in general the quantity p in equation (1) may be a function of ¢ and 
of ¢ explicitly (Rashevsky, 1940 a, 1940 b). If d¢/dt is given experi- 
mentally we may determine from equation (4), by means of equation 
(3), the functional relation between p and ¢ or t. 
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The author is indebted to Mr. H. D. Landahl for checking the 
calculations. 
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SOME CONSIDERATIONS ON MATHEMATICAL 
MOLECULAR BIOPHYSICS 


N. RASHEVSKY 
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The possibility of a single molecule or of a few molecules con- 
trolling the basic metabolic reactions in a cell is discussed from the 
molecular-kinetic point of view. For such a control by a single mole- 
cule to be possible, it is necessary to assume a chain reaction, con- 
sisting of at least two steps. 


Hitherto the development of mathematical biophysics has been 
based almost exclusively on “macroscopic” concepts of classical phys- 
ics. Diffusion phenomena, plastic deformations, and other phenomena 
important for mathematical biophysics, are treated from the point of 
view of a continuum. Different metabolic reactions are considered as 
given, and although definite assumptions are made and studied about 
the interrelations of such reactions, nevertheless these assumptions 
are of such a nature as not to involve any molecular or atomic con- 
siderations. 

It is very significant that even with this limitation to the domain 
of the classical physics of the continuum, mathematical biophysics has 
made considerable progress in quantitatively describing a large array 
of biological phenomena. Nevertheless it has been clear from the very 
beginning of mathematical biophysics, that eventually an extension 
into the domain of molecular or atomic physics will have to be made 
(Rashevsky, 1934). Even long before the development of mathemat- 
ical biophysics proper, a number of biological observations have been 
pointing to the importance for biology of atomic physics, and have 
turned the attention of noted quantum physicists to biological phe- 
nomena (Jordan, 1939). Observations on the killing of unicellular 
organisms by different radiations indicate that sometimes a single 
X-ray quantum, or a single a-particle hitting a cell are capable of 
causing its death. Although for a biologist even a large cell is still a 
“microscopic” object, from the point of view of molecular physics an 
average cell is decidedly “matter in bulk”, and at first sight it might 
look puzzling as to how a single molecule coming from outside can 
totally upset a mechanism composed of some 10% other molecules, a 
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mechanism known to follow in many respects the laws of the physics 
of the continuum. ; 

Another set of biological facts, falling into a similar category 1s 
supplied by observations on the extreme sensitivity of some sense or- 
gans, particularly that of smell. The number of molecules leaving 
some solid aromatic substances is very small in itself. If we consider 
the amazing sensitivity of the sense of smell of some animals, we may 
wonder whether only a few molecules may not be responsible here for 
setting into operation a complex macroscopic mechanism of sensa- 
tion. 

A still different field of biology leading us to the same questions 
is the study of some vitamins and hormones. Even now in many cases 
biological tests are the only ones that are sensitive enough to reveal 
the presence of those substances, the most refined physicochemical 
methods still failing. The more or less tacit assumption is frequently 
made, that surface action may account for the strong activities of 
small quantities of such substances, a very small volume concentra- 
tion being sufficient to form monomolecular surface layers, which af- 
fect the surface properties of cells. In some cases however, a simple 
calculation shows that the concentrations of the active substances are 
so small, as to be quite insufficient to cover any appreciable area of a 
cell surface. 

Phenomena of immunization, sensitization and allergy also raise 
questions of the same kind. 

It is therefore indicated to investigate first in its general aspects 
the problem of how such effects may be brought about, and to ask 
ourselves whether such phenomena do present any particular diffi- 
culties from the point of view of molecular kinetics. 

The problem consists of two parts. First there is the question 
of the atomic mechanism of catalytic interactions. To this field some 
contributions have already been made by quantum mechanics, and im- 
portant further progress is undoubtedly still to come. Second — there 
is the purely kinetic problem, as to the conditions permitting the 
above mentioned atomic mechanism to come into action. In order for 
a molecule A to affect catalytically other molecules B,, B, --- ete., 
it is necessary that a relatively close contact between the molecule A 
and the molecules B,, B, --- ete. should be established for at least a 
short time. All atomic forces, decrease very rapidly with increasing 
distance, and any distant actions are out of question. 

Consider, for the sake of definiteness, a reaction in which a mole- 
cule B breaks up into molecules C and D. Let the rate of breakdown 
of B be gq gm-cm= sec. Denoting the molecular weight of B by M, 
Avogadro’s number by N and the volume of the cell by V, we find that 
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qNV/M (1) 


molecules per cell per second are broken up. Suppose that molecules 
A of a catalyst are present in the cell, catalyzing the breakdown re- 
action. Then each molecule of the catalyst must collide every second 
with at least 


qNV/Mn (2) 


molecules B . 
Let the average concentration of B in the cell be ¢ gm-cm*, or 
éN/M molecules per cm*. Considering dilute solutions, in which the 


~ molecules behave like those of an ideal gas, denoting by o, and oz the 


radii of the molecules A and B respectively, and by v1 and v2 their 
respective average thermal velocities and putting 


__ 01 T Oe 


Ore 9D ’ (3) 


we have approximately (Jaeger, 1926) for the number Z of collisions 
which one molecule A makes with molecules B, the expression: 


_ 8u.2 + 0,2 néNo? 


Byis M ed) 


Z 


Expression (4) holds for v2 > v1, which is equivalent to the assump- 
tion that the molecule A is heavier than B. This is the more plausible 
case. When v. < v,, the quantities v, and v, have to be interchanged 
in equation (4). (Jaeger, 1926). 

The presence of the solvent will affect the number Z, but will not 
change its order of magnitude much (Jaeger, 1926). 

In order that the n molecules A could maintain the necessary rate 
of reaction q, it is necessary that Z should be at least equal to or 
larger than the expression (2). In other words, putting 


2) 2 
Stat Tous aT (5) 
BV 
we must have 
unten = qv C (6) 
- But (Rashevsky, 1940, p. 21) 
6 AQ; (7) 


ws tae ie ak Oa 


‘where ¢, is the external concentration of the metabolite B, and A is a 


quantity determined by the size of the cell, its permeability to B and 
the diffusion coefficient of B in the cell. Equations (6) and (7) give 
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Medi Rich (8) 


Ne — 
man qd 


Considering metabolic reactions with a high rate, such as the 
respiratory reactions, we have —q ~ 10-° gm cm® sec™. The quantity 
v is of the order of 10* cm sec* ; c, is about 10° gm cm, and for an 
average cell V ~ 1077 cm*. shane «© 10-%cm, and putting n=1, we 
find from expression (8): 

A 2 10‘ sec. (9) 


Such a value of A, although physically quite possible, is absurd 
for such values of g as we have assumed. It would make ¢ negative, 
since q < 0. 

If we consider the similar case of a substance B produced, then 
such a value of A will mean é~ 10-? gm cm", which is also too high. 
The values of A for different respiratory metabolites are of the order 
of magnitude of a few seconds (Rashevsky 1940). In order to satisfy 
relation (8) under those.conditions, we should have n © 10%. This still 
represents an exceedingly low concentration of the catalyst, corres- 
ponding, for molecular weights of the order of one hundred, to 10-*° 
gm-cm“. A single molecule however could not control the whole 
reaction in this case. 

Consider however a more complicated scheme. Let the reaction 
considered here, and going at a rate of 10-°° gm cm® sec be controlled 
by a catalyst, present in a concentration ¢,=—10- gmcm-, sufficient 
for the purpose. Let, however, this concentration é, be determined by 
the material balance between the rate of production of the catalyst A, 
and its rate of diffusion from the cell. Let the rate q. of production of 
A. be very small, and its external concentration be zero. In that case 
the condition that the production of A, would be controlled by a single 
molecule A would again require that A., the value of A for the sub- 
stance A,, should be of the order of 104 sec. This would require a very 
low permeability and diffusion coefficient for A.. If g-= 107° gm cm= 
sec., this will just give the correct value for é, . 

It is thus possible for'a single molecule in a cell to control a 
reaction going on at a rate of 10° gm cm sec, and involving con- 
centrations of 10° gm cm-*. But this control must be made through at 
least one intermediate stage. The single molecule A controls the 
formation of A., which goes on at a very small rate (107° gm cm? 
sec’) and is present in small concentrations (10-° gm cm-). The sub- 
stance A, in its turn controls the rate of reaction of B. When the 

“master molecule” A is destroyed, the reaction B stops due to the dis- 
appearance of A.. The concentration of the latter however will de- 
crease exponentially, with time, and the time lag between the destruc-. 
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tion of the molecule A and the cessation of reaction B will be of the 
order of A,., that is of the order of hours (Rashevsky 1940). This 
time-lag will vary as the volume of the cell. For small bacteria with 
V~10™ cm, that time lag may be reduced to only a few seconds. 
The assumption of a rather large size o, for the master molecule A, 
will also reduce A,. 

The author is indebted to Mr. H. D. Landahl for a discussion and 
check of calculations. 
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ERRATA 


In the paper “Nerve Conduction Theory: 
Some Mathematical Consequences of Bernstein’s 
Model” by F. Offner, A. Weinberg, and G. Young, 
pp. 89-103, replace s and s) by 1/s and 1/8), re- 
spectively, on page 98, line 32 and on page 99, 
line 21; replace V by V’ on page 101, line 7. Equa- 
tion (24), page 99 should read 


v=kyf(d) X d/F (d) 
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